Singular Points of Ordinary Differential Equations 

Lecture Eight 

Irregular Singular Points of Ordinary 
Differential Equations 

Solutions expanded around an irregular singular point are distinctive in one aspect: they are 
usually in the form of an exponential function times a Frobenius series. Due to the factor of the 
exponential function, a solution near an irregular singular point behaves very differently from that 
near a regular singular point. It may blow up exponentially, or vanish exponentially, or oscillate 
wildly. 

Let us start with the discussion of irregular singular points with the simpler case of the first-order 
hnear and homogeneous equation 

y' + p{x)y = 0. (6.40) 

As we know, the solution of (6.40) is 
where 



y{x) = ce-^W^ (6.42) 



P(x) = j p(x)dx 



with c a constant. If p{x) has a pole of order (A: + 1) at x = and hence has the Laurent series 
expansion 



p{x) = + ■■■ + +ao + aix^ + 



find the solution of (6.40). 
Answer 
We have 



P{x) =-^ + ---+bi]nx. 



Thus the solution of (6.40) is of the form 

exp(^-^ + ••• jx-^iM(x), 

where M(x) is a Maclaurin series. Note that x"^'M(x) is a Frobenius series. 

If p{x) has a pole of order (/: + 1) at xo instead of at the origin, and if we are interested in the 
behavior of the solution near xq, then we make the change of variable 

X - x-xo. 

The solution near xo is given by the one given above with x replaced by X. 
We are now ready to discuss the second-order Unear homogeneous differential equations (6.38) 
with an irregular singular point. 

Letxo be an irregular singular point of (6.38). If 

(x - xo)*^^c(x) = Co + ci(x - Xo) + • • • 

and 

(x - xo)^'^'^^<i(x) = do + di(x- Xq) + • • •, 

are both convergent Taylor series, and if at least one of co and do are not zero, then xo is called an 
irregular singular point of rank k. If the order of the pole of c(x) and that of d{x) at xo are {ki + \) 
and {2k2 + 2) respectively, with ki not equal to k2, then the rank k is equal to the greater of ki and 

k2. 

Note that if k = 0, then xo is a regular singular point of (6.38). 

We mention that, if xo is an irregular singular point of rank k, where A is a positive integer, then 
the solutions of (6.38) are of the form 

y{x) = exp[F(x)]F(x), (6.45) 

where 
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pf \ _ ^k-i A 1 

^""^ " {x-xoV ix-xo^-' ^"^^ 

and 

Y{x) = ^a„{x-xo)"^\ao * 0,a_i = a_2 = ••• = (6.46) 

is a Frobenius series. 

As a general remark, if we wish to find the behavior of the solution of (6.38) at very large values 
of X, we make the change of variable 

^ t ■ 

When X is very large, t is very small. Indeed, x = oo corresponds to t = 0. By making an expansion of 
the solution around r = in the way we just outlined, we obtain an expansion around x = oo. 
Keeping a few terms of this expansion is often a good approximation of the solution at very large 
values of x. 

The change of the independent variable from x to ris done with 

d ^ d ^ p. d 
dx d{\) dt ' 

dx^ dt\ dt J df dt 

-y = 0. (6.54) 



and 



Thus (6.38) becomes 



d'y , ^t-c{\) dy , d{\) 
dt^ t^ dt ^ t' ^ 



We call X = 00 an ordinary point, a regular singular point, or an irregular singular point of rank k of 
eq. (6.38), respectively, if ? = is an ordinary point, a regular singular point, or an irregular singular 
point of rank k of (6.54), respectively. In particular, we see from (6.54) that the infinity is an 
irregular singular point of rank k if either d{x) blows up like x^*"^ or c(x) blows up like x*"^ as x goes 
to infinity. 

If the infinity is an irregular singular point of rank k, then the solution of (6.38) for very large x is 
of the form 

y{x) = exp[AfcX* + Ak-ix''-^ + ■■■ + Aix]F(x), (6.55) 

where 



F(x) = ^a„x " \ao * 0,a-i = a-2 = ••• =0, (6.56) 



which is easily deduced from (6.45) and (6.46) by identifying (x - xo) with t, or 1/x. Note that the 
terms in the series (6.56) are in decreasing powers of x, and hence the series is a useful 
approximation for x large. 

Problem for the Reader: 
Find the series solution of 

y<'+X-^y = Q (6.57) 

which is useful for very small values of x. 
Answer 

The point x = is an irregular singular point with rank 1. Thus the series solution near the origin 
is of the form 

y = exp(4^)n^), (6.58) 

with 

i'W = 2^a„x«+^ao * 0,a_i = a_2 = ••• = 0. (6.59) 
We substitute (6.58) into (6.57). We get 

^ _ 2A^4I1 ^ ( 2A^ ^ Ai + l^j. ^ ^ ^^^^^ 
dx X dx x^ 
The value of A 1 is determined by requiring the most divergent term in the coefficient of Y in (6.60) 
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to vanish. There is, in eq. (6.60), a term in the coefficient of Y which is a fourth-order pole at the 
point of expansion xo = 0. Thus we require 

+ 1 = 0, 

or 

A I = ±i. 

Let us choose the root 

Ai = i. 

Then (6.60) becomes 

^ dx^ x^ dx ^ x^ 
One of the solutions can be obtained by inspection. We get 

Y{x) = x. 

This result can also be obtained if we go through the grind to obtain the recurrence formula 
(n + s - l)(n + 5 - 2)a„-i = 2i{n + s - l)a„. 
From the recurrence formula with n = 0, we get 
s= I. 

Thus the recurrence formula is 
n(n - l)c„-i = linun. 

Setting n = 1, we obtain from the recurrence formula that 

ai = 0, 

and hence a2 = 03 = • • • = 0. Therefore, the series solution for Y{x) terminates after one term. 
Thus we have obtained the closed form expression of one of the independent solutions of (6.57) as 

yi(x) = xexp(^-^y 

Taking the complex conjugate yi(x), we get the second independent solution of (6.57) as 

y2(x) = xexp(-^). 

Problem for the Reader: 
Find the series solution of the Bessel equation 

y" + x-y + (1 -p^lx^)y = (6.61) 
which is useful for very large values of x. 
Answer 

We may read off the rank at infinity from the Bessel equation (6.61) directly, without changing 
the independent variable from x to t. 

As X goes to infinity, the coefficient of y goes to a constant, or x'^. Therefore, according to this 

coefficient, 
2k-2^0, 
or 

k=\. 

The coefficient of y' is x~^. Thus, according to this coefficient, 

k-\=-\, 

or 

k = Q. 

The rank of the Bessel equation aixo = <x> is the greater of the two, and is unity. 
Thus we put 

y{x) = exp(Aix)y(x) 

where 

Y{x) = ^a„x"""*,ao = l,«-i = «-2 = ••• =0. 
The Bessel equation becomes 

{D'^ + 2AiD + x-^D + A\ + \+Ailx-p^lx^)Y=Q. 
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We require the sum of the most divergent terms in the coefficient of Fin the equation above vanish. 
Asx -y CO, the coefficient of Fin the equation above is (Af + 1). Thus we get 

A I = ±i. 

This determines the two roots of A i . 
Take the root 

Ai = i. 

Then the equation for Y becomes 

F" + i2i + x-^)Y' + ii/x-pyx^)Y = 0. 

We have 

{D^+x-^D -pyx^)Y{x) = + s){n + s+l)-{n + s) - p^]a„x-"-'-^, 

and 

i2iD + ilx)Y(x) = ^[-2i{n + s) + i]a„x-"-'-\ 

Thus we have 

2i{n + s - l/2)a„ = [{n + s - 1)^ -/7^]a„-i. 

Setting n = 0, we get 

5 = 1/2. 

For n * 0, the recurrence formula is 

a„ = ^^^±i^-i^^p^-i^a„_. (6.62) 

Thus we have 

^ 1 r(n+p + l/2)r(n-p+l/2) 
" (2/)" nirO? + l/2)r(-p + 1/2) 
Therefore, one of the series solution for the Bessel equation for very large values of x is 

yiCx) - X expC^x) 2^ ^^.^^^ ^ l/2)r(-/. + l/2)n! " ^ 

The other series solution for the Bessel equation for very large values of x is 

00 

y2{x) - X expHx) 2^ (_2/;,)« r(p + l/2)r(-p + l/2)«! " ^^"^^^ 

«=0 

Equations (6.63) and (6.64) give the solutions of the Bessel equation useful when | x | is large. 
This is to be compared with the series solution (6.32) useful when | x | is small. Since Jp{x) is a 
solution of the Bessel equation, it must be a hnear superposition of yi and y2. Thus we have 

Jp(x) = cy](x) + c*j2(x). 

The coefficient c cannot be determined with the present analysis, but we will be able to show in 
Chapter 8 that 

^-i7T{l/4+p/2) 



42n 

Thus 



Jp{x) « cos(x - 7r/4 -pn/2), x » 1. 

The series (6.63) and (6.64) differ from the series (6.32) in one important aspect: the former are 
divergent for all values of x! This can be seen from (6.62), which shows that, when n is very large, 
the ratio anla„-\ goes to infinity. 

That these series are divergent does not imply that they are not useful for the purpose of 
approximation. Indeed, they are called asymptotic series and give very good approximations to the 
solutions when x is large if the number of terms is chosen appropriately. 

The classic example of an asymptotic series is 

We have 
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7(0) = e-*dt = 1. 
J 

Thus we expect 

Kx) « 1, 

when X is small. 

To seek an approximation better than this, we approximate the factor (1 + tx)~^ in the integrand 
of 7(0 by 

-r^ — « l-tx + fx^ + ■■■ + (-tx)". 
Using this approximation for the integrand of 7(x), we get 

ft QC 

I{x) « e-'[l -tx + fx^ + ■■■ + i-tx)"]dt 
J 

= 1 + + ••• +n!(-x)". (6.66) 

The first term of this series is unity, an approximation we have already obtained. However, we note 
that the series in (6.66) diverges for all values of x as we let n ->■ co. 

Nevertheless, (6.66) is a useful formula for approximation, as we shall prove. We have 

-r-^ = 1 - tx + fx^ + . . . +i-txr + ii^^)!!^, 
l+tx ^ ^ l+tx 

which is exact. Thus we have 

I(x) = e-'[l -tx + t^x^ + ■■■ + {-txY]dt + R„, (6.67) 
J 

where 

Jo l+tx 

Equation (6.67) is exact. This is to be compared with (6.66), which is obtained by dropping R„. 
How big is the term dropped? We have 

poo 

\ Rn \<\ e-'{txy^^dt = x"^\n+ 1)1 
J 

Let a; = 0.1. Then the sum of the first five terms of (6.66), which is of the order of unity, differs from 
the exact value by 7?4, which is estimated to be about one tenth of one percent. We may improve this 
approximation somewhat by taking more terms. But as the series eventually diverges, it is obviously 
not a great idea to be too enthused and take a large number of terms in (6.66). Where do we stop? 
Very roughly, we should stop when the upper bound of 7?„ no longer goes down as we increase n by 
unity. For X = 0.1, this happens when n = 10. 

If X = 0.01, and if we take one hundred terms in the series (6.57), the error is as small as lO"'^^. 



Homework Problems due Wed next week ( Oct 13, 04) 
Chapter 6, problems. 
Chapter 7, problem 1 and 2. 
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